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O ■ Abstract 

We have developed a consistent theory of the Heisenberg quantum antiferromagnet in 
the disordered phase with a short range antiferromagnetic order on the basis of the path 
integral for spin coherent states. We have presented the Lagrangian of the theory in a 
form which is explicitly invariant under rotations and found natural variables in terms of 
^ '. which one can construct a natural perturbation theory. The short wave spin fluctuations 

are similar to those in the spin wave theory and they are of the order of the parameter 
l/2s where s is the spin magnitude. The long wave spin fluctuations are governed by the 
nonlinear sigma model and are of the order of the the parameter 1/N, where N is the 
number of field components. We also have shown that the short wave spin fluctuations 
must be evaluated accurately and the continuum limit in time of the path integral must be 
performed after all summation over the frequencies u. In the framework of our approach 
P^ | we have obtained the response function for the spin fluctuations for all region of the 

frequency u and the wave vector k and have calculated the free energy of the system. We 
have also reproduced the known results for the spin correlation length in the lowest order 
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The theory of the two-dimensional Heisenberg antiferromagnet (AF) has attracted great 
interest during the last years in connection with the problem of AF fluctuations in the 
copper oxides |], Q, f|, |3j . We especially call the attention of the reader to the review |J in 
which a general situation of the quantum AF (QAF) has been elucidated. The approach 
of these papers was based on the sigma model, which describes the long wave fluctuations 
of the Heisenberg AF in the paramagnetic phase with a short range antiferromagnetic 
order ||. The sigma model is the continuous model for the unit vector n(t, r), n 2 = 1 
in the 1 + 2 time and space dimensions J7|, §]. As a long wave theory the sigma model 
can make a lot of physical predictions such as the structure of the long wave fluctuations 
and the magnitude of the correlation length 0, f|, [J. The basic parameters of the sigma 
model, the spin rigidity p s , the velocity of sound c s , and the perpendicular susceptibility 
X± = Ps/c 2 s are calculated independently in the framework of the spin wave theory ||. 

The derivation of the sigma model presented in the papers Q and ||, without taking 
into account fluctuations, is equivalent to the statement that the classical equations of 
motion for the spin fluctuations in the QAF are described by the sigma model. But up to 



now a consistent theory of the spin fluctuations for the QAF with short range AF order 
was absent. Such is precisely the topic of this paper. 

Our approach to the description of the QAF is based on the functional integral for 
the generalized partition function (GPF) in terms of spin coherent states. This approach 
permits to solve some problems of the theory. 

This paper contains the following basic results: 

1) A method of constructing spin coherent states (Appendix |A]) invariant under ro- 
tations was proposed. It permits to write out the Berry phase and its generalization for 



the final time step the Perelomov phase [12] for the QAF in a form which is explicitly 



invariant under rotations. By the way the Lagrangian for QAF is explicitly invariant 
under rotations which permits to construct the theory of QAF with short range AF order 
in an invariant manner (section |2.2| ). 

2) Variables f2 and M, which describes AF and ferromagnetic spin fluctuations, have 
been introduced . These fields ft and M are not free. Inded ft 2 = 1 and fi ■ M = 0. 
We remove the first constraint with the help of the well known method of the Lagrange 
multiplier A. For the elimination of the second constraint we have used the Faddev - 
Popov trick which previously was used in quantum electrodynamics (section |2.3| ). 

3) We have formulated he basic l/2s approximation (section |2.4fl in the framework of 
which we have calculated the basic objects of the theory: a) the Green functions of the f2, 
M, and A fields; b) the correlation length in a zero approximation; c) the spin correlation 
function as function of the frequency uo and the wave vector k (section [2.5| ) 

4) In the leading approximation in l/2s we have integrated the action over the M 
field and we have obtained (section [2.6|) some sort of the quantum lattice rotator model 
(QLRM) [|TJ, 0. We shall call this model the spin-rotator model. It can be useful in the 
calculations in the basic l/2s approximation. 

5) It was demonstrated that the continuous approach in time can not be applied for 
calculations of the corrections to the basic approximation. The free energy also can not be 
calculated on the basis of the continuous approach in time. The method of calculations at 
the final time step was developed (section |3TT| ) . After that we calculated the free energy 
(section |3T2]), and the first order corrections to the basic approximation (section |3~3h: the 
effective action, the magnon dispersion law, and the correlation length. 

6) We have performed the separation of fluctuations for the ft field by introducing 
some separation scale h/C, < A<C hn/a in momentum space. This separation of scales 
was performed on the basis of the Pauli - Villars transformation. As a result, we have 
obtain the long wave nonlinear sigma model with additional contributions from short wave 
fluctuations. We presented an arguments showing that these short wave contributions can 
be included in the renormalization of the spin stiffness and the velocity of sound. 

2 Basic continuous approximation 



2.1 Magnetic fluctuations and l/2s approximation 

We consider the spin system which is described by the following Heisenberg Hamiltonian: 

Hh^= J - E S r S r , S t ■ S, = s(s + 1), (1) 

^ i,i'=<i> 

where S; are the spin operators; the index I runs over a two-dimensional square lattice; the 
index I' runs over the nearest neighbors of the site I; J > is the exchange constant which, 
since it is positive, corresponds to the AF spin interaction; and s is the magnitude of spin. 
The most efficient method of dealing with a spin system is based on the representation 
of the GPF Z or the generating functional of the spin Green functions in the form of the 
functional integral over spin coherent states z, z* [JlO, |ll| or over the unit vector n, n 2 = 1, 
on a sphere JO: 



Z = Tr 



exp 



(-/?#)], (3 = 1/T, (2) 



oc 



Z= ■■•/ Dfi(n)exp(A(n)), (3) 

J — -oo J — oo 

D^(n) = l[ 2 -^5(nl l -l)dn Tl (4) 

tI 

where T is the temperature, r is the imaginary time, and A(n) is the action of the system. 
In the continuum approximation, which is valid in leading order over l/2s the expression 
of the action A(n)is simplified 



A(n) = - rj2 C tot(r,l)dr, C tot (r,l) = C kin (r,l) +H(r,l), (5) 

Jo i 

Js 2 
Ckin{rJ)=is(l-cos9 T i)(p T i, H(t,1) = — E n ri ■ iM', (6) 

where 8,(p are the Euler angles of the unit vector n = (cos cp sin 9, sin (p sin 9, cos 9), and 
(p T i is the time derivative of this angle. The kinetic part of the action Chin is highly 
nonlinear and it is not clear how to proceed with it consistently. Some essential steps in 
this direction where made in |1|, [J but they do not have a final character. 

Further in this paper we use the idea of the near AF order. Following this fundamental 
hypothesis, we split our square lattice into two AF sublattices a and b. In the sublattice 
a the spins S are directed along some axis f2, in the sublattice b they are directed in the 
opposite direction. In this way, we obtain a new square lattice with two spins a and b in 
the elementary cell with a volume 2a 2 , where a is the space distance between spins. The 
axes of this new lattice are rotated by 45 degrees with respect to the primary axes. We 
assume that this AF order is only defined locally and any global AF order is absent. As 
a result, the summation over the lattice sites I and I' can be expressed as a summation 
over I G a and /' G b. which specify the space positions of the spins in the sublattices a 
and b. In terms of this notation, the Lagrangian C kin can be expressed as a sum of two 
such Lagrangians, one for the sublattice a and another for the sublattice b in terms of two 
vectors n a (r,l) and n&(r, I'). The Hamiltonian 7i conserves its form if / G a and V G b 
but J/2 — > J because the double summation is absent now. 



In this way we have two spins in each AF elementary cell which are defined in the 
different space positions I and I'. This circumstance is not convenient for subsequent 
nonlinear changes of variables. One can introduce new variables n a ^(r, I) which are both 
defined at the sublattice a (or at the center of the AF elementary cell). 

For that we pass to the Fourier image n ab (r, k) of the original vectors n a (r, I) and 
n{,(r, /'), where the momentum vector k runs over the AF Brillouin band. We can return 
to the space representation and consider the coordinate p as continuous variable. As a 
result we have the following definition: 



n a ,b(r, p) = y/l/N 8 J2 ex P (* k " P) n a,b(r, k), (7) 

k 

where 2N S is the total number of sites in the space lattice. Of course, we assume periodic 
boundary conditions. We can put the variable p now on the sublattice a (or in a center 
of the AF elementary cell). One can check that the Lagrangian C^ n will be the same in 
terms of the new variables n a ^(r, p) = n &(T, I). By the same manner one can change 
the measure of integration (|J) and write out it in terms of n a tb (r, I). The Hamiltonian "H 
preserves its simple form in the momentum representation. 

2.2 Invariant Lagrangian 

The form of the Lagrangian Ckin is not invariant under rotations although the physical 
problem itself is invariant. The reason for such situation was explained in detail in the 
Appendix (|A.2|) . In the Appendix ( A.3 ) an invariant form of the Lagrangian Ckin was 



proposed (|103|) . For two sublattices a and b it has a form 

Ckin = ~isn a ■ [m a x m a ] - isn b ■ [m 6 x m 6 ]. (8) 

Eq.(||) for Ckin is valid for any choice the unit vectors m atb if they satisfy conditions: 
m 2 ,^ = 1, n a> b ■ m a ^ = 0. For the problem of the quantum AF with two sublattices we 
can choose the following expression for vectors m a ^ 

n b - (nft • n a )n a n a - (n a ■ n 6 )n 6 

m a = =-, m b = - ^=. (9) 



(n a ■ n fe ) 2 y'l - (n a • n 6 ) 2 

This choice means physically that the vectors n a , m a determines the reference frame for 
the sublattice a, and the vectors rif,, m b determines the reference frame for the sublattice 
b (see Appendix ( |A.3| )). Substituting these expressions for m ajb into Eq. (H) we get for 
Ckin an invariant over rotation form for Cki n and also Ji (see (^|)) 

C-kin = z (n ar ; — h brl ) ■ [n ar i x n^j, (10) 

1 _ n arl ■ n brl 

H = Js 2 J2 n ari ■ n M ', n aTl E a, n bTl E b. 

Now we can introduce new more convenient variables £1(t, I) and M(t, I) which realize 
the stereographic mapping of a sphere: 

±n (i - m 2 /4) -rnxMi „, 

no,6= l i + M 2/4 ' " =1 ' n - M = - ( n ) 



In terms of these variables the total Lagrangian Cq,m = C-kin + 'H has the final form 

Ckm= ITr^' ^ = ^ 2 f E > {^-^[(1-M 2 /4)(1-M' 2 /4) (12) 

-M • M'] + ft ■ M' Q' ■ M}(1 + M 2 /4) -1 (l + M'74)- 1 , 

where ft = Q T i, ft' = Cl T i>, M = M T ;, M' = M T //. After this change of variables the 
measure of integration Dji{n) (see (£§)) passes into 

D H<P) = II (2 ^ 2 1)2 ^(nL " m<i - l)dn arl dn M = (13) 



where the product in (13) is performed over the AF (doubled) lattice cells. 



2.3 Gauge transformation 

The variable f2 is responsible for the AF fluctuations and the variable M for the fer- 
romagnetic ones. The ferromagnetic fluctuations are small according to the parameter 
l/2s and therefore one can expand the Lagrangian C^,m (|12| ) over M. The vector of the 
ferromagnetic fluctuations M plays the role (up to factor 2s) of the canonical momentum 
conjugated to the canonical coordinate Q. The term of first order in M coincides (after 
change of variables) with previous results Manousakis |]] and Sachdev || . 

From Eq. (|T2]) one can easily extract the quadratic part in the variables f2 and M of 
the total lagrangian, C qua d, 

Cquad = 2is(M • SI) + J S 2 Yl [^ ~ ^ ■ n ' + M 2 + M ■ M '] > ( 14 ) 

i'e<i> 



The Lagrangian C qua d ( |T^ ) is very simple but the measure Djj, (|T3|) is not simple due to 
the presence of two delta- functions. Therefore we can not simply perform the Gaussian 
integration over the field f2 and M. 

To solve this problem we shall use the method of the Lagrange multiplier together 
with the saddle point approximation [0, ^ to eliminate 5(Q 2 — 1): 

/oo 
exp(A x )d\, (15) 

-oo 

-A x = £ x = (t\ + c,)(tl 2 -l) 

where A is the Lagrange multiplier, and c^ is a constant which will be fixed with the help 
of the saddle point condition [0, [J . 

To eliminate 5(Q ■ M) we shall use some kind of the Faddev-Popov trick which was 
proposed in [|TJ]. Let us consider the integral: 

1=1 f(M)5(ft ■ M)dM, (16) 



and insert it in the right hand side of the identity: 



(det(B gaug )) 



1/2 



1 



exp ( --ipB gaug <p ) dip/V27r, 



(17) 



where B gaug is a positive number or a positive definite operator for some multi dimensional 
generalization. After changing the order of integration over M and <p we can make the 
change of the variable M: M — > M — flip. After that, due to the delta-function, we have 
ip = (fl ■ M) and the delta-function disappears from the integral ([16]): 



/ = (det(B gaug )Y? 2 / /(M tr ) exp (A gaug ) dM/V2n 

1 



-A 



n'liifi 



£ 



gaug 



fl-M)B gaug (fl-M) 



M^ = M - fl(fl ■ M) 

With the help of the identity ([Tj|) we can remove the delta-function 5(fl ■ M) from the 
measure (|l"3D. As a result we must substitute M — > M tr in the Lagrangian Cq,m (|12D and 



add the gauge fixing Lagrangian C gaU g due to the additional exponent in (|18D . It is very 
convenient to chose the Lagrangian C gaU g i n the form 



£gaug = Js 2 £ [(fi-M) 2 + (fi-M)(fi'-M') 



(19) 



«'£</> 



Such choice kills the major dependence on fl in the Lagrangian (14) which appears due 
to substitution M — > M 4r . We can also substitute M ir — > M in the first term of the 
Lagrangian C qua d dH) due to the identity (fl • fl) — 0. In this way, the expression ([14]) 
for C qU ad is valid in the leading order with respect to l/2s. 
The final expression for the GPF Z of the Heisenberg AF is 



Dfi(fl, M, A) exp [A(fl, M, A)] , 



(20) 



Dfi(Q, M, A) = Z gaug Dfi(fl, M)Dfj,(X), Z : 



gaug 



det(B, 



gaug j 



1/2 



Dfi(fl, M) = J] (^+Jffl— ^) rff2( r imi{ r I ) 



dim{\) = n 



r/ (2tt) 3 (1 + M 2 /4) 3 
2d\(r,l) 



L T f(2 S +l)(2vr)V2' 

where Dfi(fl, M, A) is the measure of an integration, and the action A(fl, M, A) is deter- 
mined by the total Lagrangian C tot = Ccim + £gaug + £a ( |T^ , |r5| , |T9"| ) . 



2.4 Properties of the basic approximation 

One can pass to the q = (a>, k) momentum representation (u is the frequency, k is the 
wave vector), and write out the total quadratic part of the Lagrangian in the matrix form 



C quad (q) = sX* q A(q)X q , X; 



q 
fl q 



(n:,M!), 



(21) 



P' 



P k + c 



/-<• 



rf>/2J, 



(Qk, -Pk) = J (l ±7k), 7k = (l/2)(cos(k x a) +cos(k y a)), J = Jsz. 



Here X g is a two component vector field which combines the vector fields Q q and M g ; the 
constant c^ (|T5|) is expressed through the constant /io which is the mass of the f2 field in 
the lowest order of perturbation theory. One can invert the 2x2 matrix A(q) and get the 
bare Green function G q of the Q q and M g fields 



L q = u 2 + u; 2 k , uj 2 0k = P^Q k = (1 - 7 ^ 2 + (1 + TkVoA 

Here c<j ok is the primary magnon frequency in the paramagnetic phase, G^, G^, G™, G^ 1 
are notations for the matrix elements of the matrix Green function G q . 

At first let us discuss the parameter of the perturbation theory. One can see from 
an explicit form of the Lagrangian (|12|) that the spin wave nonlinearity of the theory is 
caused by the term M tr and its modifications. Its average value is 



, w _(N -1)T P k 

2S uj=27rnTM ^ + W 0k 



<M t 2 r >=(iV-i)EGf = ^— L J±- E -rhr, (23) 



where N = 3 is the number of components of the f2 field. Summation over uj is obtained 
by standard methods || and we have 

Here, iV = 3 is the number of components of fields M and ft, n ok = (exp(u;ok/^) — 1) _ 
is the Plank function, and summation over k means the normalized integration over the 
AF Brillouin band: 

£/(k)=/^/(k). (25) 

The constants Cmo and Cmoo are defined by the relations 



Ca/o = E \V~^ = c ° " c i = 0-65075, ^ = ^-^ = 1.48491, (26) 
k V 1 + 7k k l + 7k 

where all sum of the type of (^) are calculated by the following method 

E /(7k) = / /(e)p(e)de, c n = E -7=, (27) 

p(e) = E^-7k) = ^(l-6 2 ), 
k 7r 

where K(x) is the complete elliptic integral of the first kind. By the same manner one 
can calculate the different time and space average 

<M l (r + 5,\ + r)n 3 (r,\)>=5 l ~ E ^^ ~T^r = (28) 

ZS u=2irnT,k U + W 0k 

i_ ^ exp(ik-r- \S\u^) 

"4* ¥lJ t l-exp(-^ok) ' 



Here e(5) is the sign function: e(5) = 1 for 5 > 0, e(5) = —1 for 5 < 0. For r = and 
5 J < 1 we have an explicit expression 

<M i( r + M)%(r,l)>=£%w{ (r/ ^ K r>J, ' < 29 > 

From Eqs. (^4|,[28|) we clearly see that the parameter of perturbation theory is l/2s at 
low temperatures T <C J and T /{2s J) at high temperatures T > J . Thus, perturbation 
theory is working when l/2s is a small parameter and the temperature is not high. We 
remind the reader that this is just the applicability condition of the spin wave theory. 

Now we can consider the saddle point condition for the A field < Q. >= 1 which is 
the most important constraint of the theory which determines its phase state: 

q ZS LU=2iTnT k U +U; 0k 

l=<tt 2 >=^£^(l + 2n 0k ), (31) 

AS . ZCUok 



The right hand side of Eq. ( pT|) contains two terms. The first term Qk/2cc>ok is responsible 
for the quantum fluctuations of the ft field. The second term Qk^ok/^ok is responsible for 
the classical thermal fluctuations of the f2 field. The role of these two terms is quite 
different. The quantum fluctuations are small with respect to the parameter of pertur- 
bation theory l/2s and for basic approximation they can be neglected. The quantum 
fluctuation can be considered in the continuum approximation for which Qk — 2J7" and 
k-'ok — cJLk 2 + ^ 2 i where cq s is the primary velocity of sound and Cq s = 2 J 2 s 2 za 2 /h 2 . The 
integration over the two dimensional momentum k can be easily performed 0]). The 
integration over the angle is trivial, and integration over the modulus |k| is performed if 
we introduce new variable of integration x = /3wok- As a result we have the constraint 
condition 

TN TN 

.ln(l-exp(-/i/T))~--— ln(/VT)). (32) 



2ttJs 2 v " v r/ " 2irJs< 

The coefficient before the logarithm in this equation is always small when the regime of the 
weak coupling is valid. At small temperatures T <^ J this is obvious. At the temperature 
T ~ J this coefficient coincides with the parameter of perturbation theory (123) and also 



must be considered as small. This means the logarithm in (^) must be negative and large 
in modulus. This leads to the condition fi/T ^C 1 which justifies the last simplification in 
(|32|) . As a result we have the well known p], |^, |4]] zero order expression for /i 

/i = Texp( — ^H 5 $, = hc s /n, (33) 



where £ is the correlation length. From Eq. (p3|), the very important conclusion follows: 
in the regime of the weak coupling the correlation length £ is much larger than the lat- 
tice constant a. This conclusion makes possible the scale separation for the problem of 
disordered QAF §. 

8 



To close the theory it is helpful to define the polarization operator 11(g) of the f2 field 



A A9 w = -^£A*(g)n(g)A(g), 
1 q 



(34) 



which is simply a loop, and the Green function of the A field is 11(g) 1 . In the lowest 
approximation n(g) is simply a loop from two Green function G 



H ( q ) = 2NTj2G n (q')G n (q- q '). 



(35) 



Using the Green function G n (q) from (|22|) we can perform the summation over uj and 
have the expression for the simple loop 



As 2 Xf uj iU j 2 



UJ 2 + UJ± 



+ 



UJ 2 — UJ\ 



UJ 2 + (o> 2 + UJ\) 2 UJ 2 + (uj 2 - Wl 



(36) 



where the index % :— 1, 2 corresponds the momentum l^ = k' and k 2 = k — k'. The main 
contribution in l/2s in (|36|) is from the thermal fluctuations even at low temperatures T, 
because the integral strength of such fluctuations is fixed by the saddle point condition 



(30) and does not depend on temperature. The explicit form for n (g) may be obtained 
in two limiting cases hq ^> T and hq <C T, where q 2 = uj 2 + c 2 k 2 . In the first case the 
momentum k' ~ T/c s ^C q, and we can separate the integration over k' and put k' = in 
all places in (|36"D affected by the factor m = n(k'). The result is extremely simple 



(37) 



n (g)=4^(g)= 2 f ( 2 1 + 7 2 k) g»fc T , k T = T/c s 



s(uj 2 + u 2 kJ 



At small g <^ c s /a a similar result was obtained in ||. Notice, that it exceeds the quantum 
contribution in ( pB]) n (g) = A^/4g in the large parameter lQsJ'/Nq. The second limiting 
case hq <^ T lies in the purely continuum region. It corresponds to pure classical two 
dimensional case: u = 0, and u/ = 0. Integration over k' can be easily performed, the 
result coincides with [ffl up to the normalization factor: 



Ilo(fc) 



8NT 



irs 2 q\ / q 2 + Ajj? 



In 



q + Vq 2 + 4/i 2 



2/x 



Celt. 



(3? 



2.5 The spin correlation functions 



The approach of this paper allows us to find the spin correlation functions in all values of 
uj and k. The dynamical spin susceptibility is determined by the relation 



X (oj, k)<% = -- j Q drJ2Tr{ [S^r, 1), S,-(0, 0) 
■ exp w(F — H)j\ exp (iujr — zak ■ 1) , 



(39) 



where F = — Tln(Z) is the free energy, and the wave vector k runs without limitations 
over the Brillouin band. It is well known that the dynamical spin susceptibility x( w ; k) 



coincides with the temperature Green function continued on the imaginary frequency uj. 



It can be calculated on the basis of the functional integral (|20"D 

is 
x(w, tySijSiu - u')5(k - k') = -— < n*(n, M, w', k') nj (fl, M, w, k) >, (40) 

Here, the unit vector n(f2, M, u, k) is a function of the fields f2(r, /), M(r, /) according 
to (Tfl|); the brackets < ... > mean averaging over the fi(r, /), M(r, /) fields according 
to (Iffi) . Eq. ( |4~C|) reduces the problem of calculation of the spin Green function to the 
problem of the calculation of the averages of the Q and M fields. In the lowest order in 
l/2s it is sufficient to use the lowest order relation 

n(n(r, 0, M(r, I), r, I) ~ e ial ^n(r, I) - [Q(r, I) x M(r, I)], (41) 

where q^ = (ft/a, 7r/a) is the AF vector (|TT|). Substituting the vector n from pip 
into (|40|) we get the dynamical spin susceptibility as a sum of two terms xi u > k) = 
Xa(^, k)+x^(u;, k). The spin susceptibility Xa(u, k) is responsible for the AF fluctuations. 
It is proportional to the Green function G^} analytically continued and shifted by the AF 
vector q^F 

l l ^ Js 2 z{l + lk *) 

XA[v,k) = -— , 42 

2 (a; 2 — c<jg k » + iuo) 

where k* = k— <\af, ^ok * s the rnagnon frequency ([22]). For the spin susceptibility Xa(u, k) 
we have a loop expression 

Xf(w, k) = 2s 2 £ [^(g')G M (g - g') - G«(g')G d (g - g')] . (43) 

In the main approximation in l/2s the first term gives the main contribution in the same 
manner as in (J35|) when the contribution from the Plank function m was dominated 



XF(^,k)~-^n!-^-^- 



UJ 2 +UJ\ 0J2 — Wl 



a; 2 — (u; 2 + ^i) 2 + iw^ ^" 2 — (u>2 — oji) 2 + zcj5 



(44) 



where the notation is the same as in (|3l)|). In the case of g = wcj 2 — c 2 A; 2 ^> kx the 
expression ( (44|) is substantially simplified on the basis of the idea of dominant small 
k' ~ k T , 

, ,n 2s 2 ^,,. . Js 2 z(l-7 k ) ._. 

iv A/(cu z — Co>o k + zcuo) 

For the case q < kx the expression for xf{w, k) is not so simple and we present the result 
the limit /i <C u;, Wk "C &t 

a 2 T 2 a 2 T 

Xf (u;, k) = 5-^(2) + ie(u) ■ (46) 

Trcjq 167rc;|g| 2 

■ (27r% 2 V(4^ - 3|g|) + e(-q 2 )T\q\Z(2)) , 

where £(n) is the Riemann zeta function, and e(u) is the sign function. The two terms 
in fl46| ) which are proportional to £(2) are generated by the integral which is cut at large 
k by the Plank distribution function n(k'). The ferromagnetic spin susceptibility Q47D is 
suppressed in comparison with the antiferromagnetic one (f42j) by the parameter (qa) 2 . 

10 



2.6 The spin-rotator model 

In the basic approximation in l/2s the Lagrangian is quadratic in the M field and one 
can integrate over the field M and obtain the final action for the field f2 and the Lagrange 
multiplier A: 

£nx = s(fiQ _1 n) + s(QP'Q) + +i\ (ft 2 - l) , (47) 



where the quantities Q and P' are defined in the k representation in fl22|) . As a result 
of such integration the M field is a function of the f2 field: M = Q~ X Q,. One can easily 
recognize in ( ^3|) the Lagrangian of some kind of the Quantum Lattice Rotator Model 
(QLRM)||. However it is different from the standard model due to the momentum 
dependence of the kinetic term in Eq. (0). We shall call such kind of models spin- 
rotator (SR) models. The SR model describes a quantum antiferromagnet in the limit 
s — > oo. The QLRM is also well defined and does not contain any divergences. It allows 
to perform all calculations accurately because all physical quantities are well defined in 
the framework of this model. 

3 Beyond the continuum approximation in time 

3.1 Basic approach 

When we try to construct the perturbation corrections to the basic approach discussed 
in the previous section we meet a fundamental difficulty: the integrals arising from the 
Green functions (22), over the frequency u, are not well defined. This is obvious from 



the consideration of the average < Mj(r + S, l)Qj(r, 1) >. The result essentially depends 
on the time shift 5 fl28| ) which reflects the phase space nature of the f2 and M variables: 
the Green functions G u,d (q) ~ uj^ 1 at large u. Moreover there are some doubts that the 
quantity < M ( 2 r > was calculated completely correctly. Actually, at low temperature 



[ 2 

l tr 



e r. ^4^# - o M0 , ( « 



-7T/A 27T UJ 2 + U), 



Ok 



where A is the time step in the accurate definition of the GPF presented in the Appendix 
A (0), the interval (— 7r/A,7r/A) gives the one dimensional Brillouin band for the final 
time step A. The limit A — > gives a correct value of this average. At first sight this limit 
is trivial because the integral over u in fl48|) is well defined. Suppose that we were not 
so accurate and the numerator of ( f4"8| ) in fact contains a small corrections J{\ — 7k) =>- 
i7"(l — 7k) + Atu 2 . If now we at first calculate the integral fl48|) with this numerator and 
after that pass to the limit A — ► the result will be different: Cmo —> Cmq + 1- This 
example should convinced the reader that for a spin system the continuum limit of the 
path integral for the GPS (@-§) must be produced with sufficient accuracy: at first it 
is necessary to formulate the theory at finite A and one can put A = only after the 
calculation of all integrals over u has been done. Notice also that the accurate version, 
discrete in time, of the GPS (0-fD is necessary when we calculate the free energy. 
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Now, on the basis of the results of the Appendix A we get an accurate expression for 
the quadratic part of the action. Instead of the expression (||) for the action A(n) we 
shall use a more accurate expression 

N T 

M") = - E E A [^in(j, I) + H(j, I)} , (49) 

j=0 I 

where r = jA, and AN T = (3. According to the results of the Appendix A, Ckin{j,l) 
consist of two parts Ckm — £>mod + C-pha- The first term is pure real the second term is 
pure imaginary. 

According to Eq. (|98|) the Lagrangian C mo d can presented for two sublattices a and b 
in the form 

A£ mod = -s In [(1 + n„ ■ n )(l + n 6 • n 6 )/4] , (50) 

where n a = n a (j, /), n a = n a (j + 1, /), n fe = n&(j + 1, I), n^ = rifc(j, /). It is assumed that 
vectors n a , n a , n b , n^ are functions of the dynamical variables f2 and M according to 
Eq. dTJ). 

The Lagrangian C p h a is not so simple and according to Appendix A (|106|) it is 

5^raJ' (51) 



where the quantity T a ^ is defined in ( |108| ) and has a rather complicated form. The 
expansion of AC p h a over the field M contains only odd powers of M. 

The Hamiltonian 7i(n) can be obtained on the basis of Eq. (|95D for the matrix element 
of the spin operator S if we substitute them in the Heisenberg Hamiltonian 

H(n) = Js 2 E S(n,n).S(n',n'), (52) 

Z'e</> 

where n = n a (j + 1, /), n = n a (j, /), n' = n fe (j + 1, /'), n' = n b (j, I'). All these vectors are 
also the functions of the dynamical variables f2 and M according to Eq. ([TT|) . The vector 
«S(n, n) is determined by the relation 

<n|S|n> n + n-i[nxn] 

<S(n,n) = : = -, (53) 

v ' <n|n> 1 + nn v ' 

and it transforms as a vector under rotations according to the discussion in Appendix A. 
Expanding the Lagrangians C mo d (0)' ^P ha dUD? anc l ^ ne Hamiltonian (|52] ) in powers 
of the vector M up to second order we get 

A£ kin = s[l - n ■ ft + M 2 - M ■ M + z(0 ■ M - ft ■ M)], (54) 

H = Js 2 E [Q ■ n ~ n ■ n ' + M • M + M • M' - z(Q • M - n ■ M)], 
i'e<i> 

where the usual quantities f2 and M are for the arguments j, l; the underlined ones are 
for j + 1, /; with prime are for j, /'; the underlined with prime are for j + 1, 1'. 
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According to the analysis performed in the section ( |2.3| ) it is necessary to add to the 
Lagrangian ([54] ) the gauge Lagrangian C gaug generalizing ([19]) in the case of finite time 
step 

ACgaug = 8 [<fl • M) 2 - (fi • M) (« ' M)] + (55) 

+AJs 2 ^ [(fl-M)(^-M) + (0-M)(ri'-M / )]. 
Z'e</> 

This Lagrangian C gaug kills the most strong interaction between the f2 and M fields. 
In this step we can pass to the q = (u;,k) representation 



X(g)) = p/N t J2 e wk ''- !Auj X(j, /), N t = N T N S (56) 

3,1 

u := 2nTj, (-JV T + l)/2 < j < (N T - l)/2, 

27rm a 27rm b \ N„ b - 1 A" aib - 1 



k= lWW>>' — T~ <^< 2 ' 

where we choose N T , N a , Nf, as odd numbers, N s = N a Nb is the number of sites on one 
sublattice, j,m a ,mb are natural numbers. 

Now one can write the action in the form A(ft, M) = J2 q £quad{<l) where the Lagrangian 
Cquadiq) is presented in ([H|) but the matrix A(g) acting on X g = (Q g , M g ) is now different 

k{q) = ( -l«fq), P(q))' AM (<?) = 1 - c - + A ^ + 7k), (57) 

A n (q) = l-c w + AJ( CuJ - 7k ) + A / i 2 /2 l 7, A u (g) = s w (l - AJ) 

where c^ = cos(clA) and s^ = sin(coA). Inverting the matrix A(q) we get the Green 



function generalizing (p2[) in the case of finite time step A 

-i 1 / A M (g), -A u (q) 



G « = i Afe) ) =SZwU*h). A"(i) J- (58) 

L(q) ~(1-AJ + A^/4J)[2(1 - c„) + A 2 u, 2 k ], 
where the quantity Qk, -Pk> an d the bare frequency cj k were defined in (p2|) . At small a; 



when cuA <C 1 if we neglect small terms of order Aj7" and A/i 2 ,/^ the matrices A(q) and 
G g pass (up to normalization factor A and A" 1 ) into their continuum analogues (|STJ) and 

Of course, the difference between the continuum expressions ( pil) and fl22"|) for A(g) 
and G g and the precise values ([57]) and (|58| ) are only essential for the intermediate steps 
of the calculations. We want to stress that the amplitudes of the magnon scattering in the 
skeleton approximation are determined by the continuum action (p0|) and the continuum 



Green function fl22|) only. 

At first we demonstrate that the simple averages calculated in the section ( |2.4j ) in 
fact were not calculated with sufficient accuracy. The tricks necessary to carry out the 
calculation are discussed in the Appendices. The result at low T <C J are presented 
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below. 



< MiMj >= 4^u(l + Co - ci), 



< M[M 3 >-. 



aS S ^ 

As 13 ' 



+ Ci - C 2 ) 



< MiMj > 



< M^ >= —Sijico - ci), 

<M i M' >=— ^-(ci-ca) 
4s 

4s 



< fijM,- > 



-» i,J 3 



9 "U' 

7T^S 



N 



iV\ 



< Q 'fij >= T7M 1 - 7"(1 - — + c o - 02)), < OA' >= T?%(! - 77 



N 



< n&j >-- 



4s 

N 



TT 



N 



As' 



Ar -^'( 1 -^( 1 + c o-c 2 )), 



< QiQj >- 



1a 

TV % '' 



(59) 



This calculation was performed with formulae similar to (2^.28.3^), but the expression 
(|58|) was used for the Green function G, and summation over u was restricted by N T 



terms ( 111|) . Notice, that the presence of the member 1 — c^ in the numerator of the 
Green function fl5"5|) leads to the difference of the underlined and not underlined averages. 
This difference is caused by the contribution of large uo ~ n/A. In particularly the average 

in fact coincides with the average < M ir • M 4r > but the actual 

|) on the constant (N — 1)/As. 



< M ( 2 r > calculated in 

value of < Mj r > is different and differs from the result of 



3.2 Free energy 

After the formulation of the theory for the finite time step A one can calculate the GPF 
and the free energy of the QAF in the paramagnetic phase. We can perform the calculation 
in the basic l/2s approximation. The free energy has three contributions as it follows 
from Eq. © for the GPF Z 



F AF = -T ln(Z) = F nM + F x + F t 



gang- 



(60) 



In the lowest approximation in l/2s, Zqm, Z\, and Z, 
The explicit form of these determinants follows from fl54j,|57 



TNN, 



CIM 



£ln£(?)], *I 



gaug 



i<jk 



7W. 



£ln[s 2 n (g)], 



i^'k 



,, n „g are powers of determinants. 
, and 1 



35|) 



TN X 



£ln[2sQ(g)], 



n;k 



(61) 



where all notation is in ([57]). Let us consider these three free energies separately. One 
can check that Fqm has finite limit at A — > 0, AN T = (5. F gaug and F\ do not have finite 
limit at A — > 0, AN T = (3 separately, but their sum has a finite limit. 
Consider at first F\. We present it in a form F\ = F\i + F\h. 



TN« 



'XI 



5>[no(<z)/IIoo(<Z)], F; 



7W, 



Xh 



£ln[s 2 lWg) 



(62) 
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where Hoo(q) = 4G°(g) is the polarization operator at large frequencies u ^> T ( ftTD with 
the Green function G n (q) from Eq. fl5"8|) . The summation over uj in the Eq. (R^) for F; 



A/ 

is convergent because the function ln[n (g)/Il oo (g)] tends to zero at large frequencies u 
as 1/uj 2 . Therefor the summation over u can be extended to infinity in the limit A — > 0. 
It is reasonable to joint the free energy F\ h with F gaug 

TN - TN - 1 

5F = F Xh + F gaug = -^J2\n[2sG%)/Q(q))} = -^£ln[Z(g)] = "A (63) 

cjk wk 

The summation over the frequencies u; for Fq,m can be performed on the basis of the 
Appendix (|C|) and we have the total free energy 

Faf = {{N - l)/N)F nM + F xl , (64) 

where Fq M and F\i are 

F nM = -NN s J/2 + NN S J2 Hk/2 + T ln[l - exp(-^ ok /T)]} , (65) 



ojk 



g (^ 2 + ^ 2 k )n (g) 
2.7(1 + 7k) 



where the polarization operator Ho(q) is defined in (^6|). In this formula we neglect a 
small contribution of the order n 2 / J 2 . 

The first term in the free energy Faf (HD represents the free energy of the ordered 
antiferromagnet, which consists of the ground state energy and the free energy of the 
magnon gas with two degenerate degrees of freedom. 

The temperature dependent part of the free energy (p5|) at small temperatures T <^ J 
is proportional to Faf ~ N S T 3 / J '. Such contribution has two origins: one from Fqm and 
another one from F\i. 

3.3 Perturbative corrections 

In this section we present the result of the calculation of corrections to the mass operators 
of the Q and M fields. The detailed analysis of such corrections is far beyond the scope of 
this paper. We restrict ourselves only the lowest order of the perturbation theory in l/2s. 
In this case these corrections can be presented as renormalization of the initial quadratic 



Lagrangian (154]). It is necessary to have the Lagrangian C mo d and the Hamiltonian Ti up 
to fourth order in the field M, and the Lagrangian C p h a up to third order. 
From the Eq. (p0|) for C mo d we have 

AC ( ll d = s [2K 2 - (M 2 ) 2 - M 2 M 2 + 2(M • M) 2 + (66) 

2K(M 2 + M 2 - 2M • M) - 2(0 • M - O • M) 2 } /8, K = I - ft ■ Q. 

From the Eq. ( |5T] ) for C pha we have 

ACfh = isU[(Q ■ M)M 2 - (O ■ M)M 2 - 2U(l - K)(Q ■ M - £l ■ M) ■ (67) 

•(M 2 + M 2 - 2(M • M))]/2 + 2isU 3 (n M - £L ■ M) 3 /8, U = (1 + ft ■ O)- 1 . 

15 




T<0 C=^ 



a 



Figure 1: Feynmann diagrams for the averaging of four Q fields. The solid lines correspond 
the Green functions G of the Q field, the dotted lines correspond to the Green functions 
of the A field. The small boxes are vertices which joint two Q fields in some manner. In 
this sense the diagram a describes the Gaussian contribution, the diagram b describes a 
correlation, the diagrams c and d describe corrections due to the interaction with the A 
field. The diagram c is a small correction to the diagram b, but the diagram d is of the 
same order of magnitude as the diagram b. For the case of averaging < fi^f^ >= 1 the 
diagram d simply cancels the diagram b. Notice, that the diagrams b and d are of the 
order 1/N with respect to the diagram a. 

The expression for the Hamiltonian Ti fl52|) depends on the scalar products of the fields 
to, to, to! , to', M, M, M', M'. The four order over M it is rather cumbersome and we 
will not present its explicit form. We only note that it is a pure algebraic problem. It 
is sufficient to substitute in Eq. (|53| ) for the matrix element of the spin operator S the 
expression ( |TTD for the vectors n ab . After that the result must be substituted into the 
expression for the Hamiltonian 7i (|52[). After expanding this expression over the field M 
up to fourth order we get 7i(4). 

On this step one can perform the averaging of the Lagrangians C mod , C p l a , and the 
Hamiltonian 7i( 4 ) over the fields to and M according to the rules (^). In this point 
we notice that averaging the fourth and higher powers of the to field is a little more 
sophisticated. It is necessary to take into account the fluctuation of the A field in the 
skeleton approximation over it as it shown on Fig. |I[ 

To avoid this complication we have presented the result in 1/N approximation where 
this complication is not essential. 

The effective kinetic Lagrangian and the Hamiltonian in the first l/2s approximation 
are 

A£ kin = s[a (l - to ■ to) + b (M 2 - M ■ M) - ie (to ■ M - to ■ M)], (68) 

n = Js 2 Y, M 1 -to-to) + o 2 (i -to-to')+ o 3 (i - to' ■ to) + hM 2 + 

i'e<i> 
+b 2 M ■ M + 6 3 M • M' + 6 4 M' ■M-ie 1 (to-M-to-M)-ie 2 (to'-M-to- M% 

where the notation is the same as in (0), and the constants oo, ..., e 2 are 

a i = a° i + ga h h = b° t + g(3 h e i = e° i + gj h (69) 

og = l, 6{j = l, e° = l, oS = -l, $=1, a° 3 = 0, 

b° 1 = 0, b° 2 = l, b° 3 = l, b° 4 = 0, e? = -l, e° = 0, 
a = 2, p = (3 - c + ci)/2, 70 = (7 - c + ci)/4, 
a x = 2(1 + d — c 2 ), a 2 = 3 + 8/n 2 - 3c + 2ci + c 2 , 
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a 3 = -4/tt 2 + 2c - 2ci - c 2 , (3 X = (5 - 12/vr 2 + c x - c 2 )/2, 

fa = (-4 + 8/vr 2 + c - 5ci + 4c 2 )/2, /3 3 = (2 - 24/vr 2 - 3c + ci + 2c 2 )/2, 

ft = 3/vr 2 , 7l = (-6 + 16/vr 2 + 9c - Cl - 8c 2 )/4, 72 = 1, 

where (7 = (N — l)/4s, i := 0, 1, 2, 3. The reason why the number of the components N 
enters in the effective coupling constant as N — 1 is as follows. The short range fluctuations 
are directed perpendicular to the long wave fluctuations and their number of independent 
components is iV — 1. 

Now one can write out the effective quadratic form for the Lagrangian QB8D in the u, k 
representation. Its matrix elements according to the notation ([57]) are 



A n {q) = (1 - Cuj )[a + AJ(a x + a 3 7k)] + ^Ja 23 (l - 7k ) + Afi 2 /2J, (70) 
A M (q) = (1 - c w )[6 - AJ(b 2 + 6 47k )] + AJ(b 12 + 6 347k ), 
A u (g) = ^[e + AJ(ei + e 27k )], A d (q) = -A u (g), 

where a^- = Oj + o^ and so on. The Green function is obviously represented by Eq. ( [58] ) 
with L(q) in a form 

L(g) = (e 2 + Aj7a; 1 )[2(l - cj) + (1 - c w ) 2 AJx 2 + A 2 .; 2 ], (71) 

^k = J 2 [e 2 a23(l - 7k) + Ai 2 /(J r2 6i 234 )](&i2 + &347k) 

where 6i 234 = b 12 + 6 34 , and the explicit form of the coefficients Xi, x 2 follows from (|70|). 
From (|TTIBB[[7D|,I7TJ) one can find the spin correlation function in this order in l/2s. For 
that it is necessary to take into account nonlinear corrections to the Eq. fl4"I|) which follows 
from ( ]TT| ) and also corrections to the Green functions in the framework of this formula. 

We shall give the explicit result for the correlation radius in this order in l/2s on the 
basis of Eq. (|30"D . The contribution of different frequencies u and momenta k in this 
constraint relation can be separated into two parts. The first part is the high frequency 
and momentum part. To calculation this contribution it is sufficient to take the Green 
function G (q) in the bare approximation ( |58| ) because this contribution is of the order 
l/2s. The second contribution which is proportional to the distribution function n-^ can be 
considered in the continuum approximation but with l/2s corrections taken into account: 

g"(g)- 9 2A 2^ 2 > X±=P./& ^ = C 2 k 2 + /, 2 , (72) 

where p s = Js 2 a 23 , and c 2 = eo" 2 a 23 6i 234l 7 2 a 2 /,z. Now, instead of the Eq. ( pT[ ) we have 

(N/Asp s )J2— = R, R=l-g(l + c + c 1 ), (73) 

k ^k 

The factor R inclyudes in itself the direct short wave renormalizations. Performing the 
integration in the same manner as in ([!2|) we have 

^ = Texp f — T7\r) ' Ps=PsR, £ = hc 9 /fi. (74) 

The actual temperature dependence is changed in the pre exponent factor (T —>■ J) if we 
take into account the long wave fluctuations in the next order in 1/N approximation M. 
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4 Separation of scales: description of long— wave fluc- 
tuations 

In some situations one can separate the long wave and short wave spin fluctuations. Our 
consideration of this problem has a qualitative character and we only sketch the possible 
approach to it. 

The fluctuations of the M field are short wave and the Green functions G M (q), G u,d {q) 
are regular in the long wave limit. Therefore the separation of scales is actual for the Q 
field or the Green function G n (q). This dressed the Green function G n (q) determines the 
action of the long-wave sigma model. This action is universal and can easily be obtained 
from the SR Lagrangian (^7[) or from the Green function fl58|) by the naive long-wave 
limit fj: 

A = J C cont dTd\ C cont = %± (ri 2 + cKdmdin) + ^V) + iX (n 2 - l) , (75) 

where di is the space derivative for i := x,y; c s is the velocity of sound: c 2 = 2J 2 S 2 za 2 ; 
X± = 1/2 Jza 2 is the transverse susceptibility; fi 2 is the mass of the n-field in a disordered 
phase 0. Because the characteristic low-energy scale max(T, //) is much less than the 
exchange constant J, the long-wave AF fluctuations contain many universal properties 
Q. The connection of the constants which determine these universal properties with 
parameters of the original Heisenberg model was obtained not by a direct manner but by 
some calculations in the ordered phase || ^j. 

Our basic idea is to separate scales in the Lagrangian C tot (|49|). The simplest idea 
is to separate all perturbation theory integrals over k into two parts by some separation 
scale A. But this method introduces some arbitrariness. We will use the scale separation 
based on some sort of Pauli-Villars transformation. The following identity holds: 

1/2 r d(f exp(-l<pG- 1 <p-V{<p))= (76) 



(det(G) 



IX, 



( 27 r)V2^V 2 

-i/2 f°° d(pidtp 2 ( 1 a_i 1 



(det(G , iG 2 )) j y 2 exp (--<fxG 1 Vi - ^G 2 l (p 2 - V(<pi + <p 2 )J , 

where G = Gi + G 2 , and V((p) is an arbitrary function such that the integrals in (|76| ) exist. 
This identity can be easily proved if we introduce for the second integral new variables 
ip = <pi + (p 2 and ip = <pi — <f 2- The integration over ip is Gaussian and can be easily 
performed. As a result we have the first integral. The identity (|76f) may be interpreted 
as follows. Let us choose G = p~ 2 and G x = (p 2 (l +p 2 /A 2 ))^ 1 , G 2 = (p 2 + A 2 ) -1 . We 
have achieved the splitting of the Green function of a scalar field ip into two parts with 
the help of the Pauli-Villars transformation. In fact this is a strict integral identity. 

Let us apply this method to the Lagrangian C tot ( f£5|) . We represent the field ft as a 
sum two new fields n and v: f2 = n + v. The operators G = G n (q), G\ = G n (q), and 
G*2 = G v (q), in this case in the (u>, k) representation may be chosen in the form (|57|,p^) : 

(G n (q))- 1 = 2sL(q)/A M (q) = F(q)L(q), (77) 

L(q) = 2(1 - cJ) + A 2 o; 2 k , F(q) = 2s/(l - Cu) + AJ(1 + 7 k)), 
(G^q))- 1 = F(q)L(q)(l + L(q)/A 2 ), (G^q))' 1 = F(q)(L(q) + A 2 ) 
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where the primary dispersion law u;ok is defined in (22). We omit in ( |77|) some small terms 



of the order A which can be essential only at calculation of the primary free energy ( |3.2|) . 
We assume that max(T, fi) <C A <C J . Only in this case the cut-off momentum A has a 
clear meaning. 

The fluctuations of the field n are long wave and the low frequency. This is motivated 
by the rapid fall of the Green function G n (q) with increase of the three dimensional 
momentum q = (uj, c s k), q 2 = u 2 + c 2 k 2 . As a result the Lagrangian Lq.\ reduces to C nv \ 
and has the form: 

Cnvx = (s/2J) (n • (q 2 + ^ 2 )(1 + <f/A 2 )n) + (78) 

(l/2)v ■ F(q)(L(q) + A 2 )v + t\(n + v) 2 , 

The long wave and the low frequency fluctuations of the field v are suppressed due to its 
big mass A, and the long wave and only low frequency fluctuations of the fields Q and A 
are essential. One can check that the long wave theory for the fields Q and A practically 
does not contain ultraviolet divergencies due to the cut off A. Actually the Green function 
of the A field is determined by the polarization operator n (g) (|37D and in the continuum 
limit has a form 

G\q) = (rio(g)r 1 = sq 2 /AJ, q > k T . (79) 

Now we the have following large momentum behavior of the elements of the diagram 
technique G n (q) ~ g~ 4 , G x (q) ~ q 2 , and T(q) ~ q°, where T(q) is the vertex. The 
only divergent diagram is the diagram for the mass operator in the lowest order but it 
is naturally subtracted one time and after that it is convergent. An additional order 
of perturbation theory leads to: (a) two additional Green functions G n (q), (b) additional 
Green function G x (q) (c) additional vertex T(q), (d) additional integration over q. As 
a result we have the renormalization factor R(q) = {G n (q)) 2 G x (q)T(q)q 3 ~ q~ 3 and the 
general convergence is improved. This means that the Pauli - Villars regularization is 
working. It is necessary to stress that without regularization factor A 2 /(q 2 + A 2 ) in the 
n field Green function we have R ~ q. This means that the original long wave and low 
frequency theory is unrenormalizable, and it is not possible to include all divergences in 
the finite number of objects of the theory. Of course, the parameter A is an artificial one 
and must be canceled when we calculate any observable properties due the compensation 
of the dependence on A from long wave and short wave contributions. 

Let us demonstrate how it is working in some basic example. Let us consider the 
most important constraint of the theory which determines its phase state: < ft 2 >= 1. 
Substituting f2 = n + v we have 

<n 2 >=l-<v 2 >, <n 2 >= NY,G n {q) = NT £ ^G n (g) (80) 

q uj=2ttjT k 

/7r/A (jij 



•s(9 ! + M 2 )(9 2 + A 2 )' x ' ' 2s[2A- 2 (l - c„) + w 2 k + A 2 ] ' 
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Here the summation over j for the Green function G n (q) is performed in the limits ±00. 
The main contributions in the integrals in (jBOj) from the Green function G n are from 
momentum k < A, and from the Green function G v are from momentum k > A. This 
property is general for all integrals with G n , G" . The left hand side of Eq. (|80D can be 
calculated if we reformulate accurately the summation over the frequencies uj = 2nj and 
the integration over the momentum k can be extended to infinity (see |1], Q). When we 
calculate the right hand side of Eq. (|80|) we can put the temperature T equal to zero 
and replace the summation over uj by an integration and the integration over uj can be 
easily performed. The integration over k can be performed treating A/ J <C 1 as a small 
parameter. In the part of the integral depending on A, the integration over k can be 
extended to infinity. The other part is independent of A and can be calculated in the 
same manner as the constant cmo in Eq. (pip. As a result, the constraint ( |8"0"D has the 
form 

r (A + 2Tln(/i/T)) = l + -^— -— (l + Co - Cl ). (81) 



AttJs 2 ■ \r-i jj AtcJs 2 As 

We can see that dependence on A is canceled in both sides of Eq. (^) and for /i we have 
the expression (|33|) if we take into account the main order contribution in l/2s. 

Eq. (^) demonstrates some general properties of the QAF with separation of scales 
(<a. The role of the classical and quantum spin fluctuations is essentially different. 
The separation of the spin fluctuation into quantum and classical ones in perturbation 
theory is determined by the summation over the frequencies uj. The Sommerfeld - Watson 
transformation gives us the characteristic contributions proportional to 1 + 2n(k), where 
n (k) is the Plank distribution function, and determines the separation of fluctuations. 
The quantum fluctuations are independent on the temperature T. From the point of 
view of the long wave theory the quantum fluctuations are divergent in the ultraviolet 
region. The main contribution into the quantum fluctuations comes from the region in 
the momentum space ka ~ 1. The contribution of the order k ~ A are small. Actually 
the second term in the right hand side of ( [ST] ) is less than the last one by the small 
parameter hc s A/ 'J ' . This illustrate the idea that the quantum fluctuations in the long 
wave region are not essential and can be neglected. In this long wave region only the 
classical fluctuations with the Plank distribution function are essential. In this way we 
arrive to the " renormalized classical region" picture @, |J for the classical fluctuations. 
The classical long wave fluctuations possesses some nice properties. 

1) They are convergent in the ultraviolet region due to a natural cut off at k ~ kr due 
to the Plank functions 

2) The parameters which determine these fluctuations are the renormalized p s , c 2 ., 
X± = Ps/c 2 s that entered in the action ( ]75| ) 

After this discussion the formula fl74|) for the mass /i of the Q field became obvious: 
all the quantum renormalizations are included in the renormalized spin stiffness p s and 
simple loop calculation similar to (|7|,[73|) gives this result. 

The interaction between the long wave n and the A fields continue to be essential 
which leads to the modification of the simple result (|73|) |2|, |4| . 
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5 Discussion 

Our approach to the theory of QAF is based on the explicitly rotationally invariant 
formulation of the pass integral in terms of the spin coherent states. A nontrivial choice 
of variables permits to formulate the theory based on the saddle point approximation for 
the field f2, M, and A. As a result we can find the spin fluctuations at short and at long 
distances, and at short and at long times. 

We work out the method of construction of the perturbation theory over this saddle 
point. It is not trivial because demands an accurate limit to zero over the time step A. 
This permits to make the theory in some sense similar to the spin wave theory for the 
calculation of the short distance and short time fluctuations, and similar to the sigma 
model when we consider the long distance and long time fluctuations. As a result we 
calculated the free energy and the first corrections to the Green functions of the theory. 

We performed the separate scales with the help of Pauli-Villars transformation. This 
permits to separate out the quantum and the thermal fluctuations. The quantum fluctu- 
ations are short waves and the thermal fluctuations are long waves. The fluctuations of 
the Lagrange multiplier A contains as the short wave as well as long wave contributions. 

We believe that the present approach to the QAF will be fruitful in the theory of the 
AF fluctuations in HTSC superconductors. 
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A Invariant coherent states for the rotational group 

A.l Basic formulae for the partition function in terms of the 
spin coherent states 

The partition function for a spin(|3|) system can be presented as a functional integral over 
the spin coherent states [|K], [11] 





(82) 


JVt— 1 

A(z) = -A J2 [C km (j,z)+n(j,z)}, 

3=0 


(83) 


( 1 + 2^4-1 Zn \ , N < Z* , 1 \H \Zj > 


(84) 


~ KinKJ, ~ J V, ! + kjl 2 >/ ' ^' <z* +1 \zj> 
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Here A is the time step in the imaginary time r = Aj, N T is the number of the time steps: 
AN T = /3\ Zj is the complex variable for the time step j; the integration is performed over 
z' = Re(z) and z" = Im(,2). It is supposed that the periodic boundary conditions are 
satisfied z^ T = z ; A(z) is the action, £(r) is the Lagrangian, TL(z) is the Hamiltonian. 
The vector \z > is the spin coherent state which is determined in terms of the reference 
state \ref >= \ss > 



\z >= (1 + \z\ ) exp(zS-)\ref >, S z \sm >= m\sm >, (85) 

The coherent state ( p5|) is determined up to the phase factor exp(i$(^)), where $(z) 
is an arbitrary function z, which has no influence on the partition function Z, but the 
kinetic term Cki n {j,z) is changed as the time component of the Abelian gauge field by 
the introduction of this phase factor 

£kin{h z) => C kin {j, z) - i${zj) + i$(z j+ i). (86) 

The complex variable z has the geometric interpretation of stereographic projection 
of a sphere of the radius 1/2 

* = taa(0/2)expM, \z? = \ Z ^ ) , (87) 

1 + cos 

1 — \z\ 2 2\z\ 
cos^ = —nr, sin^ = — - -, tan(<^) = Im(z)/Re(z), 

1 + \Zf 1 + \z\ 2 

where 9, (p are the polar and azimuthal angles which determine the unit vector n = 
(cos^sin^, sinyjsin^, cos#) which specifies the position on the sphere of radius 1/2. 

In the continuum limit Zj+i — Zj + Azj, where the point denotes the derivative with 
respect to time r, and the kinetic part of the Lagrangian is equal to 

z* z 
£kin(j,z) = 2s = -s(l - cos 6) <f. (88) 

1 + \z\ z 

Besides that, in the continuum limit the Hamiltonian ( p2[ ) is determined by the matrix 
element of the spin operator S 

<^|S|2>=sn. (89) 

^From expressions ( |38"| , |S9"D Eq.(0) for Ckin(T,l) and Ti(r,l) easily follows. 

A. 2 Transformation properties of the spin coherent states 

Further we want to clarify the reason why for the problem which possesses obvious rota- 
tional symmetry for the Heisenberg Hamiltonian ([]]) the Lagrangian (Q) is not explicitly 
rotationally invariant. The reason lays in the transformation properties of the state \z > 
under rotations. For further discussion is more convenient to use the state In > which 



coincides with the state \z > up to a phase factor [10, |Tl 



|n >= exp(isf)\z >= exp(—i(pS z ) exp(—iQS y )\ref >, (90) 
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in Eq.([K]) is supposed the complex number z and the unit vector n are connected by 
Eg. (|8?D . One can easily check that the state |n > satisfies the equation 

(S-n)|n>=s|n>. (91) 

Eg. (pOD can be considered as solution of Eg. (pl|) . Eg. (|9T| ) is explicitly rotational invariant 
but its solution is not. This circumstance was clarified by Perelomov in his paper in which 
the coherent states where introduced for an arbitrary Lee grope [T^]. The rotational grope 
has three parameters: Euler angles ip, 6, if) but point on the sphere is parameterized by 
two parameters (p, 9 only. Due to that there is an invariance up to the phase factor 

|n' >= exp(i$(n, a)U(a)\n >, for rii = a^rij. (92) 

Here a is the 3x3 orthogonal matrix of the three dimensional rotations; $(n, a) is the 
phase function which depends on the vector n and the matrix a; U(a) is the operator 
(matrix) of the 2s + 1 dimensional representation of the rotation a. Eg. ( J92|) clarifies the 
situation with an invariance over the three dimensional rotations. The vector |n > is not 
an invariant object but the projector |n >< n| is an invariant object. Now it is obvious 
why the total action (|83f) and the partition function ( J32|) are invariant over rotations. 
Notice that an invariance of the measure of the integration over z is obvious if we pass to 
variable n. But the kinetic part of the Lagrangian £kin{j, n) is not invariant over rotations 

AC kin (j, n) => ACkinti, n ) - ^( n j> o) + 2$(n j+ i, a). (93) 

The transformation (^3) is some kind of the Abelian lattice gauge field transformation. 
It is similar to Eq. (p6|) but has different meaning. 

^From this discussion it is obvious that the module of the scalar product of two states 
|n > and |n > is an invariant over rotations 

| < n|n > | = ( ^j . (94) 

Nontrivial invariant over rotation scalar, vector, tensor is determined by the ratio 

_ <n|4„„|n> 
Alh -~ <n|n> ' (J5j 

where Ay v .. is some operator acting in the spin Hilbert space with definite transformation 
properties. When the rotational transformation is performed the additional phase factor 
from numerator and denominator cancel each other. The invariant vector one can get if 
we choose instead the operator Ay,... the spin operator S, see Eg. (|53|). 

A. 3 Invariant coherent states 

It seems that a noninvariance over rotations of the spin coherent states |n > is their an 
integral part. However, it is possible to define formally invariant states if we introduce 
some additional unit vector m, m 2 = 1 which is orthogonal to the vector n, n • m = 0. 
After that it is natural to introduce the third unit vector k = [n x m]. This three unit 
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vectors determine a reference frame. One can determine the reference values of these unit 
vectors or the initial reference frame: m = (1, 0, 0), k = (0, 1, 0), n = (0, 0, 1). Instead 
of the transformation ( |90"1) which define the state |n > we define the state |n; m > with 
the help of the general rotation from the reference state 

|n;m >= exp(— i(pS z ) exp(— i6S y ). exp(— iipS z )\ref > (96) 

Due to our choice of the reference state \ref >= \ss > the last operator exponent factor 
in (|96|) in fact is a numerical factor. One can easily find that tan?/' = —k z /m z . For that 
it sufficient to find three dimensional matrix which rotates the vectors no, nio into the 
vectors n, m. 

Now one can find the Lagrangian £ki n {j, n) which is simply equal to 

A£ Hn (j, n) = -ln(< n;m|n;m >) , (97) 

where we use a general notation x = x(j + 1), x = x(j). It is follows from Eq.(|94"D that 

/l + n • n\ s 
<n;m|n;m>=( J F(n,m,n,m), (98) 

^kin *~"mod T *~^phai 

/l + n ■ n\ 
A£ mod = -sin I 1 , A£ pha = -ln(F(n,m, n, m)) , 

where |F(n, m, n, m)| = 1 is a pure phase factor. The Lagrangian £ mo d is pure real, and 
the Lagrangian £ p h a is pure imaginary. 

An explicit form of Y can be found from its in terms of the Euler angels 

Y 2 = exp(is(ip -V)) < ref\esxp(i9§ v )exp(i(<p - ip)S z )exp(-i6S y )\ref > /(H.C.). (99) 



The matrix element in (|99|) can be calculated using Wigner d-functions (T^, |T^] or with 
the help of z-representation of the spin coherent states (BSJ) , [TO, ITT 



a/2 
Y = P^P^PpPn, P v = [ J | "J ) > ( 10 °) 



Ti X Vlly 



P 



m z + ik z \ s fl + n-n + n z + n z + i\nxn] 



i> 



m 7 — ik z I \ H.C. 



The factor Y must be an invariant over rotations but this invariability is not obvious from 
the form ( |100| ). In fact Y must be a function of the scalar products n • n, n • m, m • k and 
so on. To find such kind of function we shall use a following trick. Let us calculate the 
factor Y in a special reference frame where m = (1, 0, 0), k = (0, 1, 0), n = (0, 0, 1). In 
this reference frame P^ = P v = P n = 1 and n x = nm, n^ = n-k, m z = mn, k z = kn. 
As a result we have the following final expression for the Lagrangian L v h a 

A£ p , a = -^lnf|=-), i? = n-m + m-k. (101) 
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In ( |101| ) and all other places we suppose that a quantity x with double underline corre- 
sponds to a quantity x without underline. Another form for the Lagrangian £ p h a one can 
get if we pass to the reference frame where n = (1, 0, 0), k = (0, 1, 0), m = (0, 0, —1) 

A£ pha = — In ([(n • n + in • k) • (m • m + ik • m)- (102) 

(1 + n ■ n n ■ m — in- k) 2 • [H.c.]^ 1 

Now one can get an expression for the Lagrangian £k in in the continuous limit. The 
magnitude of C mo d is small over A. Using obvious decomposition x — x + Ax for x : = 
n, m, k we get from ( |102j ) (or technically more complicated from (|101|) ) the following very 
simple expression for C p h a 

Cpha = -isn •[mxifaj. (103) 

One can check that if we choose for the vector m 

Z Z O / \ 

m = — , , m =1, nm = (104) 



'l-nl 

we get for C p h a = is(l — cos 9)(p, 

where e z is the unit vector along the axis z. This result is well known for £ p h a - 

A. 4 Invariant Lagrangian for quantum antiferromagnet 

The invariant Lagrangian C mo d(n) and the Hamiltonian 7i(n) where presented in the text 
(|50|.|5^) but the form of the Lagrangian £ p /i a (n) is more complicated and will be discussed 
in details below. 

The situation with the Lagrangian C pha is more complicated because it is not expressed 
in terms of the vectors n only ( |101| ) but in terms of the vectors m and k also. These 
vectors can be presented in terms of the vectors f2, M in a form 

n aib = ±c m fl + s m \, m a)b = c m \ =F s m Q, k aj6 = ±[Q x 1], (105) 

l-M 2 /4 M 2 2 

Cm ~ l + M 2 /4' $m ~ l + M 2 /4' ° m + $m ~ ' 
1 = -M 4 [fix M], |1| = 1, f = s m l, g=[flxf]. 

These expressions for m and k are not analytic over the vector M: they do not have 
definite limit independent on a direction of vector M when it turns to zero. But auxiliary 
vectors f and g have definite limit at M turns to zero. Now we shall prove that C p h a 
summing up over two sublattices a and b in fact is analytic over the vector M. According 
to ( |101| ) we have 

A^. = -\ In (§f|f). (106) 
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Substituting Eqs. ( |105| ) for n, m and k in the Eq. ( |101| ) for R we get 



R a ,b = D aib + E a>b , D a>6 = c Fn [±c m (Q.l)+i(0-k)], (107) 

E a , b = -c m s m (Q ■ il) t s rn (i ■ SI) + Cm{i ■ 1) ± i(f • k) 

The quantities D ab and E ab are not analytic over M because they contains nonanalytic 
vectors 1 and k explicitly. However, one can introduce new quantities T a;b = D* ab R ab . 
Using identity D a = —Dl one can check that the expression ( |106| ) for £ p h a will not 
change if we substitute R a ,b =$■ T a>b . This property is valid due to compensation of the 
D-multipliers from sublattices a and b. Now one can check that quantities T a>b are analytic 
over M 

T a , b = A a , b + B a , b , A a , b = D* ab D a , b = i m [l - (Q • ilf - (O • f) 2 ], (108) 
B a , b = D* atb E a>b = c m {±c m c m (Q • 0)(Q • f) + c m {il ■ f)(C • f) + 

±(f • f)(n • f) ^ (o • f) ± (o ■ n)(n • f) - 
-i\cJB. • n)(Q • g) ± c m (o • f)(o • g) + c m (^ • g)]}. 



Actually, quantities A ab and S 0j {, (|108|) depend only on vectors il, f , g and quantity c m 
which all are analytic over M. 

B Calculation integrals and sum over u 

The calculation of integrals and sums over the frequency u when the fmiteness of the 
time step A is taken into consideration demands some clarification. We start discussion 
from the case of the temperature T equal zero and after that pass to the case of the finite 
temperature. 

For the case of the temperature T equal zero the sum over u over some function f(u>) 
pass into the integral over u over the one dimensional u Brillouin band inside which 
— 7r/A < ijj < 7r/A and we are interesting in the integrals 

7T/A 

f{uj)<ko/2>K. (109) 

-7T/A 

In the perturbation theory the functions f{ui) is the rational function of c^ = cos(u;A), 
i.e. it is a ration of two polynomials from c w . Such ratio obviously can be presented in a 
form 

f(uj) = P( Cuj ) + f (O, f{cj) = Q{cu)/R{c u ), (110) 

where P,Q,R are polynomials of c w , and the power of the polynomial Q is less than the 
power of the polynomial R. The polynomial P{c w ) can be easily integrated over u. The 
rational function f'(cj) possesses some important properties. 1) This function can be 
analytically continued in the complex uj plane. 2) It have some pole singularities in the 
stripe — 7r/A < 3fJ(u;) < vr/A, — oo < 3?(u;) < oo. 3) It turns to zero exponentially in upper 
-tt/A < m(w) < tt/A, < 3f(w) < oo or low -vr/A < $1(uj) < vr/A, -oo < Qf(w) < 
half stripe. 4) The function f(uS) is periodic in the u plane f(u> + 2vr/A) = f(w). 
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Im(w) 



CI 



CI 



CI 



C2 



Re(w) 

> 



C2 



C2 



Figure 2: Contours of integration in the complex u plane for the calculation of the integral 
from the function f'(c w ). Here on the picture w = uj, wb = tt/A, and ±iw fc are the typical 
poles of a function f'(c w ). The integration is performed over contour CI if the function 
f'(c w ) is decreasing at direction 5s(uj) > 0, and the integration is performed over contour 
C2 if the function f{cj) is decreasing at direction $s(uj) > 



If all these conditions are valid we can convert the integral ( |109| ) from the function 
f{cj) into the contour integral in the complex u plane. For that one can add to the 
integral ( |109[ ) the integral over two vertical half lines which are parallel to the imaginary 
uj axis. This two lines are connects points (— tt/A, ±oo) — > (— 7r/A, 0) and (vr/A, 0) —>■ 
(jr/ A, ±oo) . The choice of the half plane for these two half axes is dictated by the direction 
in which the function f(uj) fall down exponentially. Due to periodicity the function f(u>) 
the contour integrals over u along each pair of this axes compensate each other. If we 
add these two integrals for the function f(uj) to the primary integral ( |109| ), and also the 
integral over the segment (—it/ A, ±oo) — > (7r/A, ±oo) (the last one must be infinitely 
small) we shall have the closed integral in the u plane. It is equal to the sum residues of 
the function f{ui) in this half stripe (see Fig. H). 

For the case of the final temperature it is necessary to calculate the following sum 



I = T £ /(■ 



LU 



CO 



2itjT, 



n T 



{N T - l)/2, 



(HI) 



where for definiteness the number of the N T summation points is odd. Using representation 
( |11(J| ) one can easily perform the summation over u of the polynomial Pic^) and for the 
second term this sum can be converted by the standard manner into the contour integral 



r 



Ani Jc 



cot 



'(3u 



f(cj)(k), 



(112) 



where the contour of the integration is around of the segment of the real uj axis — n/ A < 
9ft (a>) < 7r/A. By the same manner as it was discussed for the case the temperature T = 
one can supplement the contour C up to two closed contours over two the half stripes 
-7r/A < 3ft(c^) < 7r/ A, 5suj > and -tt/A < 3ft(u;) < tt/A, 5suj < 0. As a result the sum 
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( |1 1 1| ) will be equal to the sum of the residues of the function 0.5 cot(0.5 / 9a;)/ / (cw) in these 
two half stripes. 

C Calculation of the free energy 

The free energy of simple oscillator is 

F = ! X> [AJ , A, = 2(1 - c) + A 2 ul (113) 

where u = 2ttjT, and — (N T — l)/2 < j < (N T — l)/2, for definiteness N T is odd. It is 
convenient to represent the eigenvalues A w in the form 

A^, = awC, fJ-w = a + Pexp(iuA), (114) 

(a,f3)~±l + Auj /2 

In terms of the quantities p, u we have for F 

F = T^ln[/iJ = TiV T ln[a]+T^ln[l + (J3/a) exp(iuA)} (115) 

The sum over u in Eq. ( |115| ) can be easily performed if we expand the logarithmic 
function in a series over x = (a/ ' (3) exp(ikjA) and take into account only the terms x p 
with p = kN T , where k is integer, which give nonzero contributions. The result for F has 
the form 

F = TN T \n[a}+T\n[l + ((3/a) NT }. (116) 



Using Eq. ( |114| ) for a, and (3 and conditions Ao; -C 1, AN T = 1/T we easily get for the 
free energy 

F = w /2 + Tln[l-exp(-a;oT)]. (117) 

It is well known expression for the free energy of the harmonic oscillator. 
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